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Abstract. 

Based on the value of the orbital eccentricity of a particle and also its proximity to 
the exact resonant orbit in a three-body system, the Pendulum Approximation (Dermott 
& Murray 1983) or the Second Fundamental Model of Resonance (Andoyer 1903; Henrard 
& Lemaitre 1983) are commonly used to study the motion of that particle near its reso- 
nance state. In this paper, we present the method of partial averaging as an analytical 
approach to study the dynamical evolution of a body near a resonance. To focus attention 
on the capabilities of this technique, a restricted, circular and planar three-body system 
is considered and the dynamics of its outer planet while captured in a resonance with the 
inner body is studied. It is shown that the first-order partially averaged system resem- 
bles a mathematical pendulum whose librational motion can be viewed as a geometrical 
interpretation of the resonance capture phenomenon. The driving force of this pendulum 
corresponds to the gravitational attraction of the inner body and its contribution, at differ- 
ent resonant states, is shown to be proportional to e s , where s is the order of the resonance 
and e is the orbital eccentricity of the outer planet. As examples of such systems, the cases 
of (1:1), (1:2) and (1:3) resonances are discussed and the results are compared with known 
planetary systems such as the Sun- Jupiter- Trojan asteroids. 
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1 INTRODUCTION 

The study of the dynamical evolution of a planetary system while captured in a res- 
onance has a long history in dynamical astronomy. Since the pioneering work of Poincare 
(1902) on the study of the near-resonance motions in a restricted three-body system by 
means of a zeroth-order resonance Hamiltonian (i.e., a Hamiltonian with no perturbing 
terms other than the resonant ones), the body of literature produced on this subject has 
become so rich and extensive that it is virtually impossible to cite all the articles here. 
Recent discoveries of extrasolar planetary systems such as Gliese 876 (Marcy et al. 2001), 
where two planets are locked in a near (2:1) commensurability, have also provided rich 
grounds for astrodynamicists to extend such studies to the systems beyond the bound- 
aries of our solar system (Laughlin & Chambers 2001; Lissauer & Rivera 2001; Murray, 
Paskowitz & Holman 2001; Snellgrove, Papaloizou & Nelson 2001; Lee & Peale 2001). 

There are two analytical approaches that are commonly taken in the study of the 
dynamics of the bodies of a three-body system while captured in a resonance. For a 
particle with high orbital eccentricity (> 0.15) in a first-order resonance, or for a particle at 
any resonance with higher orders, the Pendulum Model is used when the particle's orbit is 
sufficiently close to the real resonant location (Dermott & Murray 1983). The Hamiltonian 
Model, or as it is often called, the Second Fundamental Model of Resonance (Andoyer 
1903; Henrard & Lemaitre 1983) is usually used for small eccentricities. A comparison 
of the results of the application of these two models to the study of the motion of a test 
particle in a first-order interior resonance as well as an exterior (1 : n'),n' = 2,3,4,5 
commensurability, with the results of numerical integrations can be found in a series of 
papers by Winter & Murray (1997a&b). In these articles, as a part of their comprehensive 
project CRISS-CROSS on understanding the location and origin of chaotic regions in the 
phase space of our solar system, Winter and Murray present a detailed analysis of the 
dynamics of a test particle at resonance. 

The purpose of this paper is to present a relatively new approach, namely the method 
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of partial averaging near a resonance, to study analytically the dynamics of a system near 
a resonant state. This technique that is based on the Averaging Theorem (see Sanders 
& Verhulst 1985 and also Wiggins 1996), enables one to avoid certain complexities by 
studying the behavior of the system averaged over a fast angular variable (Melnikov 1963; 
Guckenheimer & Holmes 1983; Greenspan & Holmes 1983). It is necessary to mention 
that a complete picture of the dynamical evolution of the system can only be obtained 
by direct analysis of its equations of motion. The partially averaged system allows one to 
focus attention on the slow-changing quantities. Such an idea is commonly used in celestial 
mechanics: the Hamiltonian of the system is averaged over a fast variable and the resulting 
averaged Hamiltonian is used to study the slow dynamics of the system. A review of this 
technique can be found in the work of Ferraz-Mello (1997) and the references cited therein. 
To demonstrate the capabilities of the method of partial averaging, a restricted, circular 
and planar three-body system is considered here and the motion of its outer planet near 
an (n : n') resonance is studied. 

Although the method of partial averaging near a resonance has long been used by 
mathematicians in their studies of dynamical systems near resonances (Melnikov 1963; 
Guckenheimer & Holmes 1983; Greenspan & Holmes 1983), its application to astronom- 
ical systems is quite recent. One can find such applications in papers by this author on 
study of the dynamical evolution of a planetary system in a uniform and homogeneous 
disk of planetesimals (Haghighipour 1999, 2000) and also in a series of articles by Chicone, 
Mashhoon & Retzloff (1996a&b, 1997a&b, 1999, 2000) on their extensive study of the 
dynamics of a binary system subject to incident gravitational radiation as well as gravita- 
tional radiation damping. 

The system of interest in this paper is a hypothetical restricted, planar and circular 
three-body system consisting of a star and two planets. The equations of motion of the 
outer planet of this system are presented in section 2. Section 3 deals with the system at 
resonance. In section 4, the method of partial averaging near a resonance is applied and 
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the averaged dynamics of the outer planet, in the first order of perturbation, is studied. It 
is shown in this section that the contribution of the gravitational attraction of the inner 
planet on the averaged dynamics of the outer body at different resonant states is directly 
proportional to the orbital eccentricity of the latter with a power equal to the order of the 
resonance. As examples of such cases, the (1:1), (1:2) and (1:3) resonances are studied in 
detail and a comparison with the system of Sun- Jupiter- Trojan asteroids as an actual case 
of a near (1:1) commensurability is presented. Section 5 concludes this study by reviewing 
the results and presenting remarks on their applicability to other planetary systems. 



2 THE SYSTEM 

The system of interest in this study is a planar three-body system consisting of a 
star S and two planets Pi and P 2 where P 2 is the outer planet and Pi, the more massive 
planet, orbits S uniformly along a circular path. The orbital period of Pi is considered to 
be known and constant. It is also assumed that the mass of S is so much larger than Pi 
and Pi that the effect of their gravitational attraction on S can be neglected. 

As mentioned earlier, it is the dynamics of P 2 that is of interest here. In an inertial 
coordinate system with its origin at S and its axes on the plane of the system, the equation 
of motion of the outer planet can be written as 



m 2 -rpr + G -rrrpr r 2 + Q 7^ ^7r( r 2 - r x ) = . (1) 

at z \r 2 \ 6 \r 2 — ri\ 6 

In this equation, Q is the gravitational constant, r\ and f 2 are the position vectors of the two 
planets, mi and m 2 are their corresponding masses and M is the mass of the central star. 
For future purposes, it is more convenient to write equation (1) in a dimensionless form. 
Introducing ro and to as the quantities that carry units of length and time, respectively, 
equation (1) can be written as 
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where f , f i and t are dimensionless quantities given by r*2 = r$ r , r\ = r fi and t = tot. 
In this equation, A; = ^Mto/rjj and = m\/M. The assumption of a uniform circular 
motion for Pi allows one to set k = 1 by choosing tq = T\ and to = Ti/2tt where 7\ is the 
orbital period of P±. 

As mentioned in the previous section, we would like to analyze equation (2) using 
the method of partial averaging. To do so, it is necessary to write equation (2) in terms 
of appropriate action-angle variables. The most appropriate action-angle variables for 
this purpose are the Delaunay variables given by L = a 1 ' 2 , G = Ve = [a(l — e 2 )] 1 / 2 , £ = 
u — esmu and g = 9 — v where e = V r Ve/ sini> and a are the eccentricity and the semimajor 
axis of the osculating ellipse of P2, is its plane-polar angle, V r and Ve are its radial and 
angular momenta and v and u are its eccentric and true anomalies related as 

G 2 

r = = a (1 — e cos u) . (3) 

1 + e cos v v ' ' 

In order to write equation (2) in term of the Delaunay variables, it is more convenient to 
first write this equation in terms of V r and Vq. That is, 
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\r — ri|- 
r sin(# - 0i) , 



- cos(e - 0i)] , 



(4) 
(5) 
(6) 
(7) 



where 9\ = u}±t and ui\ = 1, is the dimensionless angular velocity of P\. In equations (4) 
to (7), the hat signs have been dropped for the sake of simplicity and the overdot indicates 



a derivative with respect to the dimensionless time i. The vector ?! in equations (6) and 
(7) is the unit vector along f\. In terms of the Delaunay variables, equations (4) to (7) 
can be written as 



G = rF e 



(8) 



L = a(l-e 2 )- 1/2 F + e (F r sinv + F e cosw)] , 



g = - [a(l - e 2 )] 



2MV2 



,_, / sin v \ ^ . 

F e — - (F r cosv - F e sinu) 

\ 1 + e cos f / 



(9) 
(10) 



1 



£ = a ~ 3/2 + -a~ 1/2 UF r cos z;-2F e sinv) + -e (F r cos 2v - F e sin 2u) - 3 F r }>,(!!) 



where 



r — cos(6' — 6*i) 



r — ri 



\r — ri \ 6 



(12) 



3 SYSTEM AT RESONANCE 

Consider an (n : n') commensurability between the angular frequency of the inner 
planet u>i, and ug = a -3 / 2 , the Keplerian frequency of the osculating ellipse of the outer 
one. That is, 

n u>i = n' uj£ . (13) 

For our restricted circular system with this resonance condition, the Keplerian frequency 
cog, and therefore, the semimajor axis of the outer planet at resonance are constant. The 
constant value of the semimajor axis, denoted by a( n:n /), corresponds to the resonant value 

1 /2 

of the action variable L as L( n:n /) = o>^ n . n ty We would like to study the dynamics of 
the outer planet when L varies in the vicinity of this value. For this purpose, we need 
to introduce an appropriate transformation that renders equations (8) to (11) in a form 
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that includes deviations of L from L( n:n ,y Let D be an action variable measuring these 
deviations. We then write, 

L = a 1 ( H ) +^D, (14) 



and 



t = ai*£)i+<P, (15) 

where ip denotes the deviations of the mean anomaly I from its Keplerian value (Appendix 
A) . Such transformations have been repeatedly used in application of Hamiltonian averag- 
ing techniques to resonant systems (Ferraz-Mello 1997). It is important to mention that 
the choice of jU 1 / 2 in equation (14) is to assure equal lowest order of perturbation for D 
and (p after writing equations (8) to (11) for the system near resonance. Details on this 
can be found in Wiggins (1996) and also in Haghighipour (2000). 

For the purpose of writing equations (8) to (11) near a resonance, it is more convenient 
to write these equations as (Appendix B) 

dH • dH 

i = -/•-«, G = -"^' (16) 

and 

h dH . dH 

where H = — \f— ri| _1 (hereafter, external Hamiltonian) represents the perturbative effect 
of the inner planet, Pi. The dynamical equations of the system near resonance can now 
be written as (Haghighipour 1999) 



» = -^w-* D jm + ^> (18) 

<p = -3^a 7 : n/) D + n \6a~ n %D^ + + O(^) , (19) 



G = -^+0(p^), (20) 
9 = ^ + 0(^). (21) 



In these equations, all partial derivatives are evaluated at (L( n:n /), G, L^ n ,^i + (p,g). The 
averaged dynamics of the system is obtained by applying the averaging technique presented 
in appendix A to equations (18) to (21). 



4 FIRST-ORDER AVERAGED SYSTEM 

As mentioned before, we would like to study the averaged dynamics of the outer 
planet near a resonance, in the first order of the perturbation parameter fi 1 / 2 . In that 
order, equations (18) to (21) are written as 

D = -» 1/2 F, ^ = -^ 2 [^], G = g = 0, (22) 

(n:n') 

where T = dH/d£, is evaluated at (L( n:n /), G, L^ n , n ,^i + (p, g). 

In this section, we apply the method of partial averaging, as described in appendix 
A, to equations (22). As mentioned in that appendix, the analysis presented there is only 
valid for systems with one angular variable. An inspection of the main dynamical equations 
of the outer planet (i.e., equations (16) and (17)) reveals that these equations along with 
#i = 1, represent the time variations of two action variables L and G and three angular 
variables £,g and 9\. As shown by equations (22), to the first order of /U 1 / 2 , g = 0. Also, 
from the resonance condition (13) and the transformation (15), the angular variables £ and 
6>i = i are related as £ = (n/n')i + (p. This relation implies that equations (22) represent 
a dynamical system with an action variable D(t) and an angular variable <p(t). These 
equations are now in the correct form for applying the partial averaging technique. Using 
formula (A7), the averaged dynamics of the outer planet, to the first order of perturbation, 
can be written as 

£ - 3 V a (n:n')^( L ( n -- n ') ' G ' ^ ' d) = > ( 23 ) 

where the overbar denotes an averaged quantity. 



To study this equation, it only remains to calculate T which requires one to express 
H in terms of the mean anomaly £. From its definition, H can be written as 



H = 



1 + r 2 - 2r cos(0 - X ) 



-1/2 



Substituting for r from equation (3) and replacing by g + v, we have 



H = - 



1 + a 2 - 2acos(£ + g - X ) 



e a 
1 



+ - ae 
4 



2acos£ + cos(2£ + sr-6>i) -3cos(#-0i) 
2a(3 - cos 11) + 4 cos(£ + g - 0i) 



- 3 cos(3£ + £ - 0i ) - cos(£ - g + 0i) 



-1/2 



+ 0(e 3 ) 



where sin v and cos u have been replaced by 

1 -e 



cos w = — e 



+ 2 (-^) J, (,e) ' 



and 



(24) 



(25) 



sin.; = (1-e 2 ) 1 / 2 ^ sin(j£) [J^tfe) - J i+1 (je)] . 
Here J is the Bessel function of order j. From equation (25), T can be written as 

T = ^(°) + e J*» + e 2 ^ 2 ) + 0(e 3 ) , 



(26) 



(27) 



where 



J*°> = a 



-3 /2 

l + a 2 -2acos(£ + 5r-0i) sin(£ + a - X 



(28) 



(29) 



a 



1 + a 2 -2acos(£ + a-0i) 
3 



-3/2 



+ 2 a 



l + a 2 -2acos(£ + 5r-0i; 



asin£ + sin(2£ + 5r - 9±) 

1 -5/2 



(30) 



sin(£ + g-6i) 2a cos £ + cos(2£ + g - 0i) - 3 cos(sr - X ) 
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and 



^ - 1" 



1 + a 2 - 2a cos(£ + g - 6>i) 



-3/2 



- 16 a 



4a sin 2£ - 4 sin(£ + # - B x ) +9 sin(3£ + g - 0i) + sin(£ - g + 0i) 

1 ~ 5 / 2 

1 + a 2 - 2a cos(£ + # - 6»i) 
20 a sin(£ + g - dx) - 14 a sin(3£ + # - 0i) 
+ 14a sin(£- a + 0i) + 16 sin2(£ + a - X ) 
- 7 sin 2(2£ + a - 0i) + 2 (7 - 4a 2 ) sin 2£ - sin 2(g - 0i) 



15 , 



1 + a 2 - 2a cos(£ + g - 6»i) 



-7/2 



(31) 



2(4a 2 + 13) sin(£ + g - 6{) + (4a 2 - 7) sin(3£ + g - 0i) 
- (4a 2 - 15) sin(£ - g + 6i) + sm(5£ + 3g - 30i) 
+ 9 sin(£ + 3g - 30 x ) + 4a sin 2(2£ + g - 6i) - 16 a sin 2£ 
-12a sin 2(^-0!) - 8 a sin2(£ + a - X ) - 6 sin3(£ + g - X ) 



To compute the averaged value of JF, we expand [1 + a 2 — 2a cos(£ + g — 6>i)] w , w 
3/2, 5/2, 7/2, using the identity 



(1 - 2£cosa + C 2 r X = ^C q A (cosa)^ ; |£| < 1 , 

9=0 

taking into account that for the outer planet a > 1. In this equation 



(32) 



n\( \ ST r(A + a)r(A + g-/i) f , 
(7g(cosa) = > 2 — cos[(g — 2/i)q;J 

fc=o a!(g-/i)! [r(A)l 



(33) 



are the Gegenbauer polynomials. To use identity (32) for the calculation of jF, one has 
to set a = £ + g — 9\ and £ = 1/a. Simplifying .F( ),.FW and J 7 ^ using equation (32), 
one will notice that these quantities will be equal to sum of terms with a general form of 
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e s cos[(q — 2h)(£ + g — 6>i)] sin[z/£ + u'{g — 6\)] where v and v' are integers. The harmonic 
nature of these terms requires that in an (n : n') resonance, in order for the formula (A7) to 
give non-zero values, (q—2h) has to be equal to one of the four integers ±(u±n') and, at the 
same time, one of the four integers ±(1/ ± n) with the same arrangement of + and - signs. 
This immediately implies that the general term e s cos[(q — 2h)(£+g — 6 1 )] sin[is£ + is'(g — 6>i)] 
will have non-zero averaged value only if 



The quantity \Au\ in the selection rule (34) is, in fact, the order of the resonance. An 
inspection of equation (28) reveals that in an exterior (n : n') resonance of order \Ai/\, 
the first fulfillment of the condition (34) by the factor sin[z/£ + v'{g — 6ij\ appears where 
s becomes equal to \ Av\. For instance, the contribution of expansion (28) to the averaged 
dynamics of the outer planet in a resonance of the form (n : n + 1) will first appear in its 
second term, JFW, and the third term of this expansion, JF( 2 ), will be the first term with a 
non-zero averaged value when the system is captured in an (n : n + 2) commensurability. 

Let us now, just as examples of the first and the second order exterior resonances, study 
the averaged system near (1:2) and (1:3) commensurabilities. In general, after replacing 
£ by its equivalent value given by equation (15) and averaging the results, the product of 
cos[((/ — 2h)(£ + g — 9i)] and sin[z/£ + v'{g — 6>i)] will produce terms that are proportional 
to sin(n / (/? + ng). For instance, for a system near a (1:2) resonance, 



A v\ = \ v — v'\ = n' — n . 



(34) 




e (l:2) 
a 2 

"(1:2) 



o-(i : 2) sin(2y? + sr) , 



(35) 



and for a (1:3) resonance, 



^(1:3) 




(1:3) [ (3/2) 

2 [ ( 1:3 ) 
J (l:3) L 



+ 3(T 



(5/2) 
(1:3) 



+ 



(1:3) 



15 




sin(3</? + g) . 



(36) 



In these equations, 
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<7(1:2) 



E 

h=0 



a: 



1:2) 



him 



i + 



2/i + 3 
h + 1 



4a (l:2) 



2/i + 5 
/i + 2 



(37) 



r (3/2) _ 
(1:3) 



= E 

h=0 



JX1 

l (l:3) 



-l 2 



9 

2 + 



'2h + 3 s 
h + 1 . 



1 + 



1 /2/i + 5 



/i + 2 



o„2 
°"(1:3) 



(38) 



r (5/2) 
(1:3) 



E 

/i=0 



r(|+j) 
<i: 3 ) w r(|) 
i 



o„2 
OU (l:3) 



7 + 



+ 



1 - 



o„2 
°"(1:3) 



/2/i + 5\ 
\h + l ) 



Ao 2 

^"(1:3) 



/* + 

2/i + 9\l (2/i + 7)(2/i + 5) 
/i + 3 (/i + 2)(/i + 1) 



(39) 



.(7/2) 
(1:3) 



= E 

h=0 



°(1:3) ^' ^(3, 



'(1:3) 



7 

4J 



+ 



+ 



1 + 15a 2 1:3) - 4flf 1:3) 



32r? 2 

0Za (l:3) 



(2/l 



3/2/i + 7 
" 4 V /i + l 

9)(2/i + 7) 



(/i + 2)(/i + l) 



(40) 



1 fir; 2 
iDa (l:3) 



o„2 
°"(1:3) 



/2/l + 13\ 

V h + 4 ) 



(2h+ ll)(2/i + 9)(2/i + 7) 
(/i + 3)(/i + 2)(/i + l) 



Because JF is proportional to sin(n / <^+n^), equation (23) can be viewed as the equation 
of a mathematical pendulum with a potential function proportional to cos(n'<p + ng) (see 
equation (A12)). Such a pendulum, with its harmonic potential, is a characteristic of 
the first-order partially averaged system near a resonance where e( n:n /) is considered to 
be constant (Sanders & Verhulst 1985; Lichtenberg & Lieberman 1992; Wiggins 1996; 
Haghighipour 1999). For instance, for the cases of (1:2) and (1:3) resonances, 
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Figure 1. Graphs of the potential function IA(\ : 2)(<{>) (top left) and its phase diagram (top right) and the 
potential function W(i : 3)(</5) (bottom left) with its associated phase diagram (bottom right) for the system 
studied by Haghighipour (1999) at (1:2) and (1:3) resonances. The scale on all vertical axes is 0.01 . The origins 
on the horizontal axes of the graphs of the (1:2) resonance have been shifted by -0.44(rad) and the corresponding- 
origins of the graphs of the (1:3) resonance have been shifted by -1.7(rad). 



W(l:2)fa) 



3e ( i :2 ) 

Z "(l:2) 



°"(i:2) cos(2</? + g) + Constant, 



(41) 



and 



W(1:3)M = 



C (l:3) 
^"(1:3) 



a 



(3/2) 
(1:3) 



, o (5/2) 
+ dC7 (l:3) + 



15 



a 



(1:3) 



r (7/2) 
(1=3) 



cos(3(/? + (7) + Constant. (42) 



Figure 1 shows the graphs of these two potential functions with their corresponding phase 
diagrams. In producing these graphs, the numerical values of the orbital eccentricity 
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Figure 2. Graphs of D(t) and (f>{t) for the systems of figure 1, partially averaged to the first order. 



and the semimajor axis of Pi have been taken from the restricted three-body system of 
Haghighipour (1999) at (1:2) and (1:3) resonances. 

The harmonic nature of the potential function of the pendulum (23) indicates that 
this pendulum can be in three dynamical states; stable equilibrium corresponding to the 
minimum of the potential function (centers C on the phase diagrams), unstable equilibrium 
corresponding to the maximum of the potential function (saddle points S) or, an oscillatory 
(librational) motion around the stable equilibrium (the orbits inside the separatrix, the 
dashed orbit that passes through the saddle point S). The resonance lock phenomenon is 
geometrically depicted by these librational motions. 

The oscillatory variations in values of (p create a harmonic behavior for the action 
variable D. From equations (22), for a system at an (n : n') resonance where f is pro- 
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portional to sm(n'(f + ng) , D will be proportional to cos(n'ip + ng). Figure 2 shows the 
graphs of D and (p against time for the systems of figure 1. As mentioned in section 3, 
D is the measure of changes in the action variable L or, in other words, an indication of 
the variations of the semimajor axis of the outer planet from its resonant value. From 
equation (A14), the width of the resonance band within which the semimajor axis of the 
outer planet varies around its resonant value is limited by the height of the separatrix and 
to the first order of perturbation can be written as 



Aa 



(n:n') 



= 4 



1/2 



(43) 



where AU( n:n ')(<p) is equal to the difference between the maximum and the minimum 
values of W (n:n /)(^). 

It is necessary to mention that the procedure presented here for expansion of [1 + a 2 — 
2acos(£ + g — 6*!)]"™ using Gegenbauer polynomials, is not valid for a (1:1) resonance. At 
this state, from equation (13), the Keplerian period of the outer planet becomes nearly 
equal to T±. That means, a(i : i) — 1 and expansion (32) is no longer applicable. However, 
it is still possible to use the method of averaging to study the dynamics of the outer planet 
near a (1:1) commensurability. The external Hamiltonian H, in this case, must be studied 
in its entirety as given by equation (25). Expanding H to the second order in eccentricity 
and integrating T in the vicinity of the (1:1) resonance, we have 



•^(1:1) =^ a (l:l) t 2 - e (l:l)] Sm <> + 9) 

~ 3a (i:i) e (i:i) sin2(<^ + g) 
45 

+ -J a (l:l) e (l:l) Sm3 (£ + 9) 



1 + a (l:l) - 2a (l:l) COS (> + #) 



1 + a (l:l) - 2a (l:l) COS^ + 9) 



i + a?!:!) -2a (1:1) cos(v? + g) 



-3/2 



-5/2 



-7/2 



(44) 
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Figure 3 Graphs of the potential function (<p) (left) and its phase diagram (right) against (p. 



The potential function associated with the first-order partially averaged system is given 
by 
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W(1:1)(¥0 



a 



(1:1) 



1 + tt (l:l) ~ 2a (l:l) «>s(p + #) 

21 



-1/2 



20(1:1) 



e (l:l) COS (^ + f) 



i + ^l:!) - 2a (l:l) COS (¥ 3 + ^) 



-3/2 



(45) 



- ^ e (l:l) Sm V + f) 



1 + a(i:i) - 2a(i:i) cos(y? + #) 



-5/2 



Figure 3 shows the graph of this potential for g = 0°. The fact that dlA^-.i)^) / d(p shows 
slight deviations from zero near the points of stable equilibrium is an indication of a slow- 
frequency librational motion in that neighborhood. To see this, let us apply this analysis 
to the system of Sun- Jupiter- Trojan asteroid. Located at the L4 and L 5 Lagrangian points 
of Jupiter's orbit, Trojan asteroids are in a near (1:1) resonance with Jupiter and have 
a librational motion with a period of approximately 148 years (Brown & Shook 1964; 
Fleming & Hamilton 2000) . Substituting for a and e in equation (45) by the values of the 
semimajor axis and the eccentricity of Trojan asteroids and Taylor expanding U(i-\){0) 
around its stable equilibrium, one will obtain a librational period of approximately 110 
years. The difference between this period and the 148 years reported in the references 
above can be attributed to several factors such as neglecting the gravitational effect of 
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Saturn, restricting Jupiter to stay on a circular orbit and also to the coupling between 
ip and the argument of the pericenter g. Studies are currently underway to extend the 
partially averaged equations of this system to the second order of perturbation where <p 
and g decouple. A closer value for the librational period above is expected in this case. 

The first-order partially averaged system near a resonance, presented by equation 
(23) is, in general, Hamiltonian. It portrays the phenomenon of capture in a resonance 
as librational motion of a pendulum. However, with respect to an arbitrary perturbation, 
this Hamiltonian system is structurally unstable. In an actual system, in order to be able 
to draw conclusions on long-term behavior of quantities such as orbital eccentricity of the 
outer planet, its angular momentum and also the precession of its orbit, it is necessary 
to extend this analysis to higher orders of the perturbation parameter \j}l 2 . Such an 
extension will allow for the time variation of the angle g to be taken into consideration. 
This will render the system of equations (18) to (21) in a set of equations with two angular 
variables, <p and g. In order to be able to apply the method of partial averaging near 
a resonance to this system, it is then necessary to introduce an averaging transformation 
that renders the equations of the system in a form that to the first order of perturbation, it 
becomes automatically equivalent to the first-order partially averaged system at resonance. 
The second-order partially averaged system is then obtained by averaging those equations 
using formula (A7). Such studies are currently in preparation for publication. 

5. SUMMARY AND CONCLUDING REMARKS 

Application of the method of partial averaging to the study of the dynamics of the 
outer body of a restricted three-body system while captured in a resonance has been 
presented here. Analysis of the first-order partially averaged system near a resonance 
has revealed that the equations of motion of the outer planet, averaged over fast periodic 
motion at resonance, resembles a mathematical pendulum. Such a pendulum analogy can 
also be found in the comprehensive study of orbital resonances among planetary satellites 
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by Peale (1986) and also in the comprehensive study of the dynamical behavior of a test 
particle near an interior as well as an exterior resonance in a restricted three-body system 
by Winter and Murray (1997a&b). 

In the analysis presented here, the driving force of the pendulum-like first-order par- 
tially averaged system is obtained from the external Hamiltonian H (equation (24)) which 
involves the gravitational effect of the inner planet on the dynamics of the outer one. In 
writing the dynamical equations of the outer planet in terms of the Delaunay variables, 
this force appears as F r and Fg in equations (8) to (11). It is important to mention that 
the form of these equations (i.e., equations (8) to (11)) is quite general and independent 
of the physical nature of the perturbation. In a system where the perturbations are non- 
Hamiltonian, extra terms will be added to the functions F r and Fg as well as equations 
(16) and (17). However, equations (8) to (11) will keep their general form (Haghighipour 
1999, 2000). The procedure presented here regarding the application of the partial aver- 
aging technique and the analysis of the pendulum-like equation are also quite general and 
can be applied to non-Hamiltonian systems in a similar fashion. In fact, one of the most 
important features of the method of partial averaging near a resonance is that it presents 
a general analytical procedure that is equally applicable to both Hamiltonian and non- 
Hamiltonian systems. As examples of the systems where the method of partial averaging 
near a resonance has been used in conjunction with non-Hamiltonian perturbations, we 
refer the reader to Chicone, Mashhoon & Retzloff (1997a&b) and Haghighipour (1999). 

As mentioned in section 4, the first-order partially averaged system near a resonance 
presents the first step in utilizing the method of averaging in the analytical study of the 
dynamics of a resonance-locked system. To this order of approximation, the argument of 
the pericenter g was assumed to be constant. In order to obtain a more comprehensive 
picture of the dynamics of a system near a resonance and the roles that perturbative effects 
play in its stability, one has to extend this analysis to higher orders of the perturbation 
parameter \i x l 2 . Such an extension is necessary to assure decoupling of the angles ip and 
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g (e.g., see Haghighipour 2000). At that stage, one can apply the averaging technique 
presented here to the dynamics of the system at the second order of perturbation by 
averaging those equations over the fast-changing angular variable cp and studying the 
pendulum-like equation of the angular variable g. Such studies are currently underway for 
the case of (1:1) resonance and their applications to Trojan asteroids. 

The choice of a restricted system as presented in this study was merely to focus 
attention on the method of partial averaging and its capabilities as another approach to 
analytical analysis of the dynamics of a system near a resonance. The analysis presented 
here is quite general and can be applied to any dynamical system at resonance. One can 
apply such analysis to the systems at interior resonances by setting £ = a in equation (32) 
and changing the selection rule (34) to 

\Au\ = n - n' . (46) 

An implication of this selection rule can be found in a recent paper by Michtchenko & 
Ferraz-Mello (2001) on analytical modeling of the Jupiter-Saturn system near their (5:2) 
resonance. They show that at the lowest order, the contribution of the resonant part of the 
disturbing function appears as the third power of eccentricity, a result that is also implied 
by the selection rule (46). 

Other interesting cases for application of the partial averaging technique are the study 
of the stability of the extrasolar planetary system Gliese 876 (Marcy et al. 2001) where its 
two planets are locked in a near (2:1) commensurability, the study of dynamical stability 
of P-type binary planetary systems where recent numerical integrations by Holman & 
Wiegert (1999) have indicated the existence of islands of instabilities for eccentric binaries 
at (n : l),n > 3, resonances and also the study of the dynamical evolution of pulsar 
planetary systems such as PSR B1257 +12 (Konacki, Maciejewski & Wolszczan 1999) and 
PSR B1620-26 (Thorsett et al. 1999; Ford et al. 2000; Ford, Kozinsky & Rasio 2000; 
Rasio 2001). It appears that gravitational radiation reaction may play a vital role in the 
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dynamical evolution of these system (Chicone, Mashhoon & Retzloff 1996 a&b, 1997 a&b, 
1999, 2000) 
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APPENDIX A: METHOD OF PARTIAL AVERAGING NEAR A RESONANCE 

A short introduction to the method of partial averaging near a resonance as used in 
this article is presented here. For more details on this technique, the reader is referred to 
Sanders & Verhulst (1985), Wiggins (1996), Haghighipour (1999, 2000) and the references 
therein. 

Partial averaging near a resonance is based on the application of the Method of Aver- 
aging to the dynamical equations of a system in the vicinity of its resonant state (Sanders 
& Verhulst 1985; Wiggins 1996). These equations are usually written in one of the several 
standard forms (Sanders & Verhulst 1985). In celestial mechanics, it is customary to write 
the dynamical equations of the system in terms of action-angle variables. 

Consider a perturbation system with an action variable B and an angular variable (3 
such that 
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B = eM(B, (3, t, e) , (Al) 

and 

P = oj (B) +eQ(B,/3, t,e) , (A2) 

where Ai and Q are periodic in time with period T and is the frequency of the unper- 
turbed (e = 0) system. At resonance T and lv are related as 

Iujq = I' uir , (A3) 

where ut is the angular frequency associated with T and / and V are positive integers. 
One can show that in the vicinity of the resonance state (A3), equations (Al) and (A2) 
can be written as (Wiggins 1996; Haghighipour 2000) 

£ = e 1 / 2 M(B ,P,t),+e£^(B ,P,t) + 0(e 3 / 2 ) (A4) 

and 



6 = e^£^(B ) + e [ Q (B ,P,t) + ^^(fib) 



+ 0(£ 3 / 2 ) , 



(A5) 



where 



£ = £- 1 / 2 (B-B ) 



Q = P-u (B )t, 



(A6) 



represent deviations of B and f3 from their resonant values Bq and uJo(Bo)t, respectively. 
Equations (A6) are, indeed, the necessary transformations for writing equations (Al) and 
(A2) in the neighborhood of the (/ : I') resonance. The averaged dynamics of the system at 
this neighborhood is obtained by averaging equations (A4) and (A5) using the averaging 
integral 

i r lT _ 

(A7) 



n = — J n[B , uj (B )t + e , t] dt 
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The first order partially averaged system near (I : V) resonance is obtained by neglect- 
ing the 0(e) terms in equations (A4) and (A5) and is given by 



E = e 1 ' 2 M(B , Q) , (A8) 

and 

& = e^S^(Bo), (A9) 

where the overbar indicates an averaged quantity. According to the Principle of Averaging 
(Saners & Verhulst 1985; Wiggins 1996; Chicone, Mashhoon & Retzloff 1997a; Haghigh- 
ipour 1999) the dynamics of the system of equations (Al) and (A2) can be approximated 
by the dynamics of the system (A8) and (A9) during the time interval e~ x l 2 t. It is neces- 
sary to emphasize that in order to be able to make such an approximation, it is required by 
the Principle of Averaging, that the main dynamical system (i.e., equations (Al) and (A2)) 
to have only one angular variable. Extension of the method of averaging to the systems 
with two or more angular variables can be found in the works of Grebenikov & Ryabov 
(1983), Arnold, Kozlov & Neishtadt (1988) and also in a recent paper by Cucu-Dumitrescu 
& Selaru (1998) on study of the equations of motion around an oblate planet. In the first 
order of perturbation, however, such an extension is not necessary. 
Introducing Ti, as 



H = e 1/2 



(AlO) 



one can show that equations (A7) and (A8) can be written as 



Equations (All) imply that TL can be considered as the Hamiltonian of the first-order 
partially averaged system at resonance. To this Hamiltonian, one can attribute a potential 

24 



function given by 



V(S) = - J M(B , Q)dS. 



(A12) 



Differentiating equations (All) with respect to t and using the Hamiltonian H, one can 
write 



e - e 



duJo 
dB 



(B ) M(B ,e) = 0. 



(A13) 



Equation (A13) can be regarded as the equation of a mathematical pendulum with Hamil- 
tonian TC and potential function V(@). The librational motion of this pendulum presents 
a geometrical interpretation for the resonance capture phenomenon. The maximum varia- 
tion of the action variable B associated with these librational motions is given by (Wiggins 
1996) 



AB = 2l2e 



dB 



1/2 

(Bo) V Max (&) - V min (&) } +0(e). (A14) 



APPENDIX B 

From definition of H, we have 



dH 
~dL 

dH 
dG 

dH 
~dJ 

dH 



1 



\r — T\\ 
1 

\r — ?i| 
1 

\r — 

1 

I — * — * I _ 

r — ri 



-0i) 


dr 
dL 


-0i) 


dr 
dG 


-0i) 


dr 
~dl 


de 




dg- 
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+ rsin(6> - 6»i) 



de_ 

dL 

m_ 

dG 



d£ 



(Bl) 

(B2) 

(B3) 
(B4) 



From these equations it is evident that one needs to compute derivatives of r and 6 with 
respect to all Delaunay variables. From equation (3), we have 



dr ( G \ 2 ( de dv 

cos v —— — e sin v 



dv \ 

dL \l + ecosvJ V dL dLJ ^ 

dr ( G \ ( G \ 2 / de dv \ 

ttt; = 2 - cos v 7—- - e sin v , (B6) 

dG Vl + ecosw/ Vl + ecosz;/ V dG dGJ v ' 

dr . du 

di =eaSmU di (B7) 

and dr/dg = 0. On the other hand, from 9 = g + v, d9/dg = 1 and the derivatives of 6 
with respect to L, G and £ will be equal to derivatives of v with respect to these variables. 
Using £ = u — e sin-u and G = L(l — e 2 ) 1 / 2 along with equation (3), the partial derivatives 
of r with respect to the Delaunay variables can be written as 



— = - a" 1/2 (2e - cos v - ecos 3 v) , (B8) 
dL e v ' 

% = ^a(l-e 2 )] 1/2 cos,, (B9) 

% = ea(l-e 2 )" 1 / 2 sinv , (BIO) 
d£ 

and the partial derivatives of 9 with respect to L, G and £ will be equal to 



d9 Gr 2 

- = eL 3 (1 _ e2) sin,(2 + ecos,), (Bll) 

1 r , ,, n i/2 rl 1 

— = --al-e 2 7 - + -7. ^ sin., B12 

dG e L Lr a(l — e 2 ) 



Replacing the derivatives of r and 9 in equations (Bl) to (B4) by their equivalent expres- 
sions given by equations (B8) to (B13), one can write 
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dH 
~dL 



r — cos(6> — 9i) 
\f — ?i| 3 



(2e — cos v — ecos 3 v) 



+ 



sin(# - 0x) 
\r — ?i| 3 



(2 + e cos v) sin v > , 



dH 1 



= ^[a(l-e 2 )] 



2MV2 



cos(# - e x 



?ll 3 



COSf 



sin(# - Q x ) 

i -* -* i ^ 
r — ri \ 6 



+ 



r a(l — e 2 ) . 



~d£ 



a(l-e 2 ) 1/2 jesini> 



r — cos(6> — 9\ 
\r — ?i| 3 



+ ^(l-e 2 ) 



sin(fl-fli) 
|r — ?i| 3 



and 



dH r sin(# - 0^ 



|r — Fil 3 



which along with equations (12) immediately result in expressions (16) and (17) 
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